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A method to predict the emergence of different kinds of ordered collective behaviors in systems of globally coupled chaotic 
maps is proposed. The method is based on the analogy between globally coupled maps and a map subjected to an external 
drive. A vector field which results from this analogy appears to govern the transient evolution of the globally coupled system. 
General forms of global couplings are considered. Some simple applications are given. 
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PACS Number(s): 05.45. +b, 02.50.-r 

Coupled map lattices constitute useful models for the 
study of spatiotemporal processes in a variety of contexts 
|jj . There has been recent interest in the investigation of 
the emergence of ordered collective behaviors (OCB) in 
systems of interacting chaotic elements by using coupled 
map lattices [|]J|^,^| . Such cooperative phenomena have 
been considered relevant in many physical and biological 
situation s P J7|,P|,P|. In particular, globally coupled maps 
(GCM) Eg] can exhibit OCB such as: a) formation of 
clusters, i.e., differentiated subsets of synchronized ele- 
ments within the network pd| ] ; b) non-statistical proper- 
ties in the fluctuations of the mean field of the ensemble 
PlP^f ; c) global quasiperiodic motion and dif- 

ferent collective phases depending on the parameters of 
the system Q . 

Much effort has been dedicated to establishing the nec- 
essary conditions for the emergence of various types of 
OCB, mostly involving direct numerical simulations on 
the whole globally coupled system. However, this direct 
procedure gives little information about the mechanism 
for the emergence of OCB. In this letter, we propose an 
alternative method which allows us to gain insight on the 
conditions for the emergence of specified types of collec- 
tive behavior. 

The procedure is suggested by the analogy between a 
globally coupled map system 

= (1 - e)f(x-) + eH(xl,x 2 t , ...,x?), (1) 

and an associated driven map 

s t+1 = (l-e)f(s t ) + eL t . (2) 

In eq. ([!]), cc™ (n = 1, . . . , N) is the state of the nth 
element at discrete time t; N is the size of the system; 
e is the coupling parameter; and f(x) describes the local 
dynamics. A commonly used form for H is the mean field 
coupling |I^,[l^,[l4| , however the global coupling function 
H t = H(xf,xf, ■ ■ ■ in eq. ([!]) is assumed to be any 
function invariant to argument permutations. In cq. (|^), 
St is the state of the driven map at discrete time t, f(st) 
is the same local dynamics than in eq. ([!]), and L t is the 
driving term. In general, L t may be any function of time. 



Here we focus on the relationship between periodically 
driven maps and the emergence of periodic and quasiperi- 
odic OCB in GCM. In this case, the asymptotic collec- 
tive behavior is characterized by periodic or quasiperiodic 
time evolution of the mean state (x) t and of the coupling 
function H t for t > t a , where t a is the transient time 
required to reach the OCB. The elements of the system 
follow a similar behavior to H t , and in general are segre- 
gated in clusters which tend to be out of phase (!(]] . 

In the case of an OCB consisting of K clusters with 
period P, the asymptotic behavior of the system can be 
characterized by a K x P matrix 



fxl 

\Xp 



K 
Xp 



(3) 



where the fcth column contains the temporal sequence 
of the P values adopted by the elements belonging to 
the fcth cluster; and the ith row displays the state of all 
clusters at time i. The time step t have been replaced 
by the index i which runs from 1 to the periodicity P. 
Therefore, the P asymptotic values $i adopted by the 
global coupling H t> t a are 





JVfc times 

where is the number of elements in the kth cluster 
and J2k^k = N. Thus, we may express the asymp- 
totic temporal behavior of H t >t a by the vector <D = 

When the driven map, eq. (|2|), is subjected to a peri- 
odic drive L t with period P, the long-term response of 
St depends on the initial condition si, and on the spe- 
cific sequence of values adopted by L t , which we will 
denote as the vector L = (L\, . . . , Li, . . . , Lp). The re- 
sponse s t (L,si) may be periodic on a region R of the 
P-dimensional space spanned by all possible vectors L. 
For a given L £ R, there are, usually, J different asymp- 
totic periodic responses of s t with period M, each one 
associated to a set of initial conditions {si}j, where 



1 



j = 1, 2, . . . , J. If we denote by a\ the asymptotic peri- 
odic response of St(L,s\) that is reached for s\ £ {si}j, 
then the possible asymptotic responses can be repre- 
sented by the J x M matrix 



(5) 



A/ 



M 



where the jth column contains the asymptotic periodic 
response of St for the initial states si £ and the ith 

row displays the state of all asymptotic responses at time 
i. As before, the time step t has been replaced by the in- 
dex i running from 1 to the periodicity M of the response. 
Obviously, this scheme is valid if all the asymptotic re- 
sponses (j = 1, . . . , J) have the same periodicity M (not 
necessarily M = P) . The region R may consist of various 
subregions R J M , where J different asymptotic responses 
of periodicity M occur. In particular, the subregion R J P 
is relevant since there the asymptotic responses of the 
driven map have the same periodicity than the drive, and 
consequently the analogy between eq (Q) and eq. (Q) can 
be established. Note that if L = the columns of the 
OCB matrix \ are contained in the response matrix <x. 
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1. Main subregions R{ { for the driven logistic map 
= 2. Parameter values are r = 1.7 and e = 0.2. 
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As an illustration, figure 1 shows the main subregions 
resulting from eq. (||) for a logistic map f(x) = 
driven by different L t with the same period 



rx 



1 

P = 2. Parameter values are r = 1.7 and e = 0.2. The 
subregions R J M indicate where J asymptotic responses of 
periodicity M occur when the components L\ and L2 of 
L lie in the range [—1,1]. The structure of this diagram is 
actually more complex; there are small zones in between 
the marked subregions where the driven map reaches a 



variety of periodic responses, (not shown in fig. 1). We 
have also observed shrimp- like structures for L\ and L2 
outside the range [—1,1], similar to those reported for 
two- parameter maps |H| . 

Now consider the evolution of the GCM towards an 
OCB of period P. During this transient time, the ele- 
ments of the system segregate in "swarms" which pro- 
gressively shrink, and eventually form clusters. Simulta- 
neously, Ht evolves towards its asymptotic periodic be- 
havior $ given by eq. (||) . The convergence of Ht towards 
$ can be represented in the same P-dimensional space 
spanned by the vectors L as a trajectory joining the point 
(Hi, .... Hp) to (Hp+i, ■ ■ ■ , H 2 p) to (H 2P+ i, . . . ,i?3p) 
... to ($!,..., $ F ) = 

Suppose that we want to find out if, for a given cou- 
pling function H t in eq.(|l|), an OCB with period P and 
K clusters can be observed. The distribution of elements 
in clusters corresponds to a partition {pi,£>2, • ■ ■ ,Pk}, 
where pt is the fraction of elements in the fcth cluster. 
Then, based on the analogies presented above, the fol- 
lowing procedure may be employed: 

a) Determine the region Rp in which the response ma- 
trix <r is K x P, by exploring the asymptotic responses 
of the associated driven map for different drives L t of the 
same period P in eq. (0). 

b) Construct the associated coupling vector = 
(Oi, . . . , Oj, . . . , Op), whose components are given by 



JVpi times 




) (6) 

n p K times 



c) In the space Rp construct the vector field 



V = L-Q. 



(7) 



If in the region Rp the vector field V is convergent 
towards the locus where V = 0, then the chosen OCB 
can take place for appropriate initial conditions in the 
GCM system. 

As shown in the examples bellow, the vector field V 
acts as an indicator of the transient trajectory of H t ; 
that is, the region of convergence of V towards V = 
is related to the basin of attraction of the chosen OCB 
in the space of the P initial values (Hi, . . . , Hp). The 
field V is particulary efficient in pointing the transient 
trajectory of Ht when the associated driven map reaches 
its asymptotic periodic response in few iterations. No- 
tice that the vector field V can be calculated without 
previous knowledge of the existence of chosen OCB in 
the corresponding GCM. Moreover, the first step of the 
procedure (and the most time consuming) must be ex- 
ecuted only once; thereafter, steps (b) and (c) can be 
performed repeatly for many different coupling functions 
H t and different partitions {pi,p 2 , ■ ■ ■ ,Pk}- 

In order to show applications of the procedure, we 
choose globally coupled logistic maps f(x) = 1 — rx 2 , 



2 



and look for an OCB consisting of two clusters of period 
two. Parameter values are fixed at r = 1.7 and e = 0.2 
in what follows. 
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FIG. 2. a) The vector field V and the trajectory of Ht for 
the arithmetic mean global coupling, b) Magnification of (a). 

As a first example, the coupling function H t is assumed 
to be the arithmetic mean; i.e., H t — jj 5Zn=i x t ~ ( x )t- 
Therefore, ©i = picrj +p 2 crf , and Vi = 0^ — Li, (i = 1,2). 
Figure 2(a) shows the vector field V{L\, L 2 ,pi) in the re- 
gion R 2 for a partition p\ = 0.477 and p 2 = 1 — p\. The 
vector field V(Li, L2,pi) is plotted as arrows of length 
proportional to \V\, direction given by tan -1 (Vi/V2), and 
origin at (L\,L 2 ). Figure 2(a) also shows the superpo- 
sition of the trajectory of the coupling function H t con- 
structed by joining (Hi,H 2 ) to (H 3 ,H^) ... to ($i, <J> 2 ). 
Random initial conditions of TV = 2000 elements uni- 



formly distributed on [— 1, 1] are used, i.e., (x)\ — 0. For 
these initial conditions, the GCM rapidly collapses into 
two clusters with the chosen partition, i.e., N\ = 954 and 
N 2 = 1046, respectively. For other initial conditions, dif- 
ferent partitions result. In those cases the vector field 
V maintains the same appearance, except that the point 
where V = moves along the dashed curve in fig. 2(a) 
as p\ varies. The two labeled values of p\ at the ends of 
the dashed curve in fig. 2(a) indicate the position of the 
convergent point V — for these critical values of p\ . An 
OCB of period two with two clusters can not be observed 
for values of pi outside this interval. Note that for the 
first few iterations, the H t transient trajectory is roughly 
suggested by the field V, since at these early stages the 
cluster are being formed. As time progresses, the vector 
field V becomes a better indicator of the evolution of the 
GCM system toward its OCB, as showed in fig. (2b). 

In a second example, we consider the geometric mean 
of the moduli of the values of the elements of the system 
as the coupling function, i.e., H t = J] 
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fore, in the region R 2 , Oi — \a. 
Figure 3 again shows the corresponding field V together 
with the trajectory of the coupling function H t . The 
same initial conditions and system size as in the first ex- 
ample have been used. In this case, the convergence to 
the OCB is faster in comparison to the first example. 



T 



Pl =0.6 



iVVY'A 

////// //V 
i i / / / / / / / / / 
i j / / / / / / / / s 

i i J J _/ / / S i/ i/ 
S s S s is 
s S s V ^ 

s is IS a~ 




/////////; ft ,^ 

//////// / f f 1 \ \ \ 

/////////// f t \ \ 

J J A <A l/i .A £ h h If i i i t I i n 



_L 



.1 .15 .2 .25 .3 .35 .4 .45 .5 

Ll . H 2 t-i 

FIG. 3. The vector field V and the trajectory of H t for the 
geometric mean of the moduli coupling. 

Finally, the method can be used to infer a global cou- 
pling function capable of producing a quasiperiodic OCB 
with two clusters. This can be achieved by modifying the 
coupling function of a known case of periodic OCB, just 
near its convergent point V = 0. Then, we may expect 
that the modified GCM will not reach its original peri- 
odic OCB, but would remain close to it. For the first 
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example, this can be done by noticing first that when L 
is near the convergent point, one of the two asymptotic 
responses of the driven map adopts a value close to 0.86 
every two time steps. Then, the coupling function of the 
first example may be modified in such a way that only 
near 0.86 the function is drastically affected. 
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FIG. 4. a) The vector field V and the trajectory of Ht 
for the modified arithmetic mean global coupling, b) The 
corresponding asymptotic return map of the mean state of 
the system. 

We have tried with coupling functions of the type: 



f 



N 



n=l 



1 — a I 1 — exp 



1 



b(x? - 0.. 



(8) 



The results for parameters values a = 1.7 and b = 10 3 in 
eq. <JsJ) are shown in figure 4(a). Additionally, the return 
map of the arithmetic mean < x >t+i vs. < x > t is 
shown in fig. 4(b). This example shows the usefulness of 
the method for designing globally coupled systems with 
specific features. Quasiperiodic OCB appears associated 
to the type of convergence of the field V occuring in this 
example. This type of convergence can occur for other 
f(x) and H t . 

A variety of clustered OCB can be predicted for cou- 
pled logistic maps with < r < 2, provided that the 
value of e is in the appropriated range, in accordance with 
fig. 3 of In the above examples, the field V can be 
represented on a plane, but for P > 3, projections of V 
can reveal its global structure. As P increases, the com- 
putation time increases potentially with P. Then, nu- 
merical methods to search for convergence towards V = 
may be used to speed the process. 

In summary, we have presented a method to study the 
emergence of ordered collective behaviors in globally cou- 
pled maps, based on the analogy of these systems with 
a driven map. A vector field defined on the space of 
the drive term in the driven map acts as an indicator of 
the evolution of the coupling function in the GCM. The 
method can be used to predict if specific types of OCB 
can take place in a GCM system and to visualize its asso- 
ciated basin of attraction. The limitation of this method 
lies in the fact that the driven map (eq. (Q)) must posses 
periodic asymptotic responses so that the matrix er can 
be defined. A large family of maps fulfill this condition. 
The examples presented here show that progress in the 
understanding of the collective behaviors of GCM can be 
made by investigating its relation with driven maps. Ex- 
tensions of this method could be applied to phenomena 
such as control of chaos, chaotic synchronization, phase 
segregations, and intermittent OCB. 
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